Subcritical modulational instability and transition to chaos from periodicity 
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This study shows that a modulationally destabilized monochromatic wave in a fluid system may 
undergo a subcritical bifurcation directly into chaos, when dissipation is weak enough. Analysis is 
made within the framework of the complex Ginzburg-Landau equations. 
PACS numbers: 05.45.-a, 05.45. Jn, 82.40.Bj 



This study concerns a subcritical transition to chaos 
directly from periodicity in a high dimensional system. 
We consider in particular a continuous physical system 
where a monochromatic wave is modulationally destabi- 
lized. The slow temporal and spatial modulation of the 
envelopes of a destabilized wave train, in a weakly non- 
linear medium, is universally described by the Ginzburg- 
Landau equations(GL). [1] Within the theoretical frame- 
work of GL, this study reports that, when the strength 
of dissipation of the medium is weak enough, the mod- 
ulated wave undergoes a subcritical transition directly 
into chaos. In the related earlier studies, it was generally 
shown that the modulationally destabilized wave would 
lead either to a supercritical transition into a 2-frequency 
motion or to a period-doubling bifurcation. [2-4] 

We start with the partial differential equation known 
as the Ginzburg-Landau equations; [1] 



* t = + (e + i)V xx - (e - i)|^| 2 * 



(1) 



, where the real positive parameter e measures the 
strength of dissipation of the system. When e = 0, the 
system becomes dissipation-free and GL is reduced to the 
well known Nonlinear Schrodinger Equation. [5] 

The GL is invariant under x — > —x, and therefore, for 
an even initial distribution under periodic boundary con- 
ditions with a period — ir/q < x < n/q, the evolution can 
be expanded in terms of cosine functions as follows, 

^(x,t) = ai(t) + 02(f) cosqx + 0.3(f) cos2ga; + • • • (2) 

, where cii(t) are complex. 

Substituting Eq.(2) into Eq.(l), we obtain, up to 02 
[7], the following set of coupled equations; [6] 
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+ (i - e)(a 1 2 a 2 * + 2|ai| 2 a 2 + -\a2\ 2 a 2 ), 



(3) 



where a'i , a'2 denote respectively the time derivative of 
the complex variables ai(f), 02(f)- 

Eq.(3) now contains two parameters e and g, which 
measure, respectively, the strength of dissipation of the 



system and the modulation wave number. When e = 0, 
the system becomes conservative with a conservative 
quantity, E — \a\\ 2 + | ct2 1 2 • This set of equations pos- 
sesses a special symmetry property, namely, they are in- 
variant under the symmetry 0.2(f) — > —02(f). 

Eq.(3) possesses the following limit cycle solution with 
frequency u> = 1; 



ai (f) 



(4) 



A linear stability analysis shows that, for small values of 
e, the periodic orbit becomes unstable for the following 
perturbation 



if 



oi(f = 0) = l,o 2 (f = 0) = 5, (5<1, 



q<qth~V2(l~e 2 ). 



(5) 



(6) 



In this study, we investigate the dynamics of the desta- 
bilized periodic orbit when dissipation is weak enough, 
say, when e = 0.1. For e = 0.1, the periodic orbit begins 
to lose its stability at q t h ~ 1-4. As q decreases from q t h 
to q = q c = 1.236, it is found that the frequency at u = 1 
starts to shift gradually from lu = 1 to lu = lu c , where the 
frequency is suddenly replaced by a broadband spectrum 
as shown in Fig. 1. Figure 2 gives the corresponding 
bifurcation diagram, where Ki, the ith maximum of the 
signal 1 02(f)), 1S plotted for all i = 1, 2, 3, .... on the verti- 
cal axis as a function of q on the horizontal axis. The re- 
gion denoted by "Pi" or "JV' in the bifurcation diagram 
corresponds to a periodic regime with a single frequency. 
The region "C" , on the other hand, is obtained from the 
motion exhibiting a broadband power spectrum. When 
the motion is in the region C, we further plotted Afj+i 
as a function of Mi for all i = 1,2, 3, where Mj now 
denotes the ith maximum of the signal 1 1 — ai (f ) | . The re- 
sult is shown in Fig. 3, where we observe a non-invertible 
one-dimensional return map, indicating that the flow is 
now in a strange attractor. Incidently, we notice that 
this 1-d map has exactly the same shape as that of the 
Lorenz attractor [8]. 

To understand the dynamics underlying the direct 
transition to chaos from periodicity, we now rewrite the 
solution as follows; 
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oi (t) = 1 - Z{t) + iW(t), 
a 2 {t)=X(t)+iY(t) (7) 

, where Z(t), W(t),X(t), Y(t) are all real. 

We then attempt to make a geometrical analysis of 
the corresponding flows in the 4-d phase space, spanned 
by X, Y, Z, and W. Since one can not visualize tra- 
jectories in such a 4-dimensional space, we further re- 
duce the dimension of the dynamics through the method 
of the Poincare's surface of section; Whenever a trajec- 
tory passes through the plane W = with a direction 
W > 0, we mark the other three coordinates (X, Y, Z) 
of the piercing point. The nth point (X n , Y n , Z n ) is then 
uniquely mapped into the next point (X n+ \, Y n+ i, Z n+ i), 
which now defines a discrete trajectory in a 3-dimensional 
space, which we refer to in this study as the Poincare 
space. 

The original limit cycle whose instability is under in- 
vestigation is then represented in this Poincare space as a 
fixed point located at the origin. Also the property that 
the system is symmetric under the reflection a 2 (t) — > 
—a 2 (t) is now converted to the reflectional symmetry of 
the Poincare space under (X, Y, Z) -> (-X, —Y, Z). The 
fixed point at the origin in the Poincare space becomes 
unstable at the threshold of linear instability q = q t h- 
Figure 4(a) shows that, when it becomes unstable, nearby 
trajectories are pushed away spirally into one of the pair 
of stable fixed points formed nearby. This type of bifur- 
cation, referred to as a symmetry breaking bifurcation, 
results from the symmetry of flows under the reflection 
(X,Y,Z) -> (-X,-Y,Z). This symmetry breaking bi- 
furcation leads to the periodic regime denoted by P 2 in 
Fig. 2. As q is further lowered to the critical value q c , 
the trajectories are all of a sudden repelled from the new 
fixed point. Figure 4(b) displays the moment of transi- 
tion, where it is seen that the neighboring trajectories of 
the new fixed point are spiraling out directly into a global 
double-wing structure. The noninvertible f-d map dis- 
played in Fig. 3, which has exactly the same shape as the 
1-d map of the Lorenz attractor, comes from this struc- 
ture. This double wing structure is therefore a strange 
attractor, which appears to have the same geometry of 
the Lorenz attractor. This provides a geometric expla- 
nation of the observed direct transition into chaos from 
periodicity. Incidcntly, readers are referred to Ref. [9] for 
a derivation of the Lorenz model for the description of an 
infinitesimally weak amplitude modulation of a periodic 
continuous wave near the threshold of the modulational 
instability. 

If it were a subcritical bifurcation, then the transition 
must also exhibit hysteresis [10]. By reversing the di- 
rection of q, we observe that the inverse transition takes 
place at q = q c ' = 1.261 instead of q c = 1.236. Such hys- 
teresis is only possible if the pair of periodic orbits and 
the chaotic attractor coexist in the range (q c , q c '). Figure 
5 actually depicts the presence of the invisible repulsive 
invariant torus (RIT) in the Poincare space, which now 
separates the two basins; one for the global chaotic at- 



tractor and the other for the stable fixed points. This 
concludes that the transition is indeed a subcritical Hopf 
bifurcation. 

Throughout the study, one may have noticed that there 
are strong geometric similarities between the discrete 
flows in the Poincare space and the continuous flows 
of the actual Lorenz system [8]. We notice, however, 
one major difference. In the Lorenz system there exists 
the homoclinic connection, which is then followed by the 
transient chaos [11] that is to be converted into a chaotic 
attractor by a crisis [12]. We found no such homoclinic 
connection leading to a transient chaos in the present 
discrete flows of the Poincare space. 

Although we have started with the continuous flow in 
a specific 4-d dynamical system, the conclusions drawn 
on the dynamics are based on the discrete flows in the re- 
duced Poincare space. The results therefore may bear a 
generic significance, and we present in Fig. 6 a schematic 
drawing of the bifurcation pattern leading to the subcrit- 
ical transition to chaos from periodicity concluded in this 
study. Actually we observed the same subcritical tran- 
sition to chaos in the continuous GL equations for the 
same initial condition. Figure 7 displays the bifurcation 
diagram found in the GL equations where the subcritical 
transition take place at q = q c = 1.31. 

To conclude, we have shown that a periodic orbit in 
a continuous system can undergo a direct transition to 
chaos when dissipation is weak enough. An analysis of 
the underlying geometry shows that the transition is ac- 
companied by the emergence of the smooth manifolds 
of the Lorenz attractor in the 3-d Poincare space. This 
study through the idea of diffcomorphisms [13], may give 
an example which helps to understand the dynamic scope 
of the Lorenz attractor in a more complex, higher dimen- 
sional system. 
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FIG. 1. Power spectra at (a) q c 
q = 1.235 < q c . 
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FIG. 4. Geometry in the 3-d Poincare mapping space. Dis- 
crete trajectories are projected on the Y — Z plane, (a) Pitch- 
fork bifurcation of the fixed point at the origin, (b) Subcrit- 
ical Hopf bifurcation. Trajectories starting near one of the 
fixed points are spiralled away into the outside global struc- 
ture. The flow in this structure produces the 1-d return map 
shown in Fig. 3. 



" ' 1.15 1.40 1-50 

q 

FIG. 2. Bifurcation diagram. Ki denotes the ith maximum 
of |a2 (t) | - The ordinate displays Ki for all i — 1,2,3, 
"Pi" denotes the periodic state with frequency lu = 1 
before destabilization. "P2" denotes a periodic state for 
<? c < q < qth ~ 1-4. At q < q c , transition to chaos("C") 
takes place. 
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FIG. 3. One-dimensional return map for q = 1.235 < q c . 
Mi denotes the ith maximum of the signal |1 — ai(t)\. The 
abscissa shows Mi, while the ordinate displays M n +i, for all 
t = l,2,3,"-. 
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FIG. 5. Hysteresis. A repulsive invariant torus(RIT) sepa- 
rates two basins; one for the stable fixed point and the other 
for the outside chaotic attractor. Trajectories starting just 
outside of the RIT are spiraling out into the global chaotic 
attractor. 
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FIG. 6. Schematic drawing for the bifurcation pattern. 
"P" denotes a periodic state and "C" a chaotic state. The 
first bifurcation (at q = qth) is a symmetry breaking bifur- 
cation. The second bifurcation(at q — q c ) is a subcritical 
bifurcation into chaos(C). Hysteresis is observed in the inter- 
val [q c ',q c ]. Following the usual convention, the dashed lines 
represents unstable states and the solid lines stable states. 
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FIG. 7. Subcritical transition to chaos from periodicity in 
the Ginzburg-Landau equations for e = 0.1. See the captions 
of Fig. 2 for Pi, P 2 and C. q c = 1, 3, q th = s/2. 



